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Abstract 

In the framework of the Tsallis nonextensive statistical mechanics we study an 
assembly of ./V spins, first in a background magnetic field, and then assuming them 
to interact via a long-range homogeneous mean field. To take into account the spin 
fluctuations the dynamical field coefficient is considered to be linearly dependent 
on the temperature. The physical quantities are evaluated using a perturbative 
expansion in the nonextensivity parameter (1 — q). The extended Curie- Weiss law 
in the mean field case has been generalized. The critical temperature and the Curie- 
Weiss constant are found to be dependent on the nonextensivity parameter (1 — q). 
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I Introduction 



A generalization of the Boltzmann-Gibbs extensive statistical mechanics was proposed by 
Tsallis [I] via a deformation of the functional form of entropy 

S = k\n q W } q G R + , (1.1) 

where k is the Boltzmann constant, and W is the weight factor. The stability of Tsallis 
entropy fll.il) is ensured by maintaining the deformation parameter q as a positive real 
number j2]. The deformed q- logarithm and its inverse the g-exponential function read 

x l ~ q — 1 i 
hi g x = - _ - , exp q (x) = [1 + (1 - q)xY~i. (1.2) 

The extensive classical Boltzmann-Gibbs entropy is recovered from (11.11) in the q — > 1 
limit. Nonextensive statistical mechanics has found applications to a wide variety of 
fields such as anomalous diffusion |5,6J, quantum information theory [7f8] , astrophysical 
problems pfTU] . and biological systems [TTfT2"] . 

The introduction of physically appropriate constraints while the nonextensive entropy 
(11. ip of the system is maximized has been achieved in two alternate ways. In the so-called 
second constraint picture [5] one works with the unnormalized g-expectation values of a 
physical variable O: 

(O)? = £ (pf^X - c( 2 >(/3) = £ (pf\p))\ (1.3) 

j j 

(2) 

where p^- \p) is the ensemble probability of the microstate j in the second constraint 
picture: 



4 (/5) 



It is well-known that in the second constraint formalism the expectation value of the unit 
operator is not preserved, and, as given in (II. 3p . it equals to the sum of the q- weights 
c (2) (/3). 

In contrast the third constraint scenario [1] employs appropriately normalized escort 
probabilities resulting in the following g-expectation values: 



where pj (/3) is the corresponding probability of the microstate j: 
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and c^ 3 \(3) is the sum of (/-weights. The thermodynamic averages of the physical quantities 
are obtained via the generalized partition function Zg 3 \(3) that relates [1] to the sum of 
the g-weights as 

Zf\p) = {t^\p))^. (1.7) 

In both the above cases the inverse temperature (3 = (kT)^ 1 is associated jl] with the 
Lagrange multiplier corresponding to the internal energy. In the light of the fact that in the 
third constraint formalism the unit operator trivially preserves its norm for an arbitrary 
q, it is considered [1] to be fully satisfactory. Fortunately, the ensemble probabilities in 
the two pictures may be interrelated jl] by the following equivalence relation: 

*?\p) = p? ) (n (i.8) 

where the general recipe for constructing the transformation jl] to the auxiliary temper- 
ature (3' reads 

(3 = (3' {P> , (1.9) 

where Uq 2 \f3) is the internal energy in the second constraint picture. Its explicit compu- 
tation pertinent to our models will be discussed later. Here we note that as a result of the 
above equivalence property the dynamical quantities obtained in the second constraint 
framework may be translated to their respective values corresponding to the choice of the 
third constraint. 

Application of nonextensive statistical mechanics to spin systems was initiated in 
|13j . where an assembly of iV noninteracting spin-| particles in a background field was 
studied in the second constraint formalism, associating the inverse temperature with the 
Lagrange multiplier corresponding to the energy. A similar study based on the third 
constraint framework was also performed [H]. The magnetic susceptibility in this model 
showed [r3]Tbi] the interesting feature referred to as dark magnetism, indicating that the 
apparent number of spins are different from the actual number of spins. Employing a 
high temperature limit, it was observed in [2] that in the domain q > 1 (q < 1) the 
effective number of spins N e g > iV (N c q < N). In the study of manganites nonextensive 
statistical mechanics has been observed [15] to fit the experimental data on magnetization 
better than the standard Boltzmann-Gibbs statistics. A numerical analysis of a multilevel 
spin model has been done [16] following Tsallis statistics. A collection of spin clusters has 
been examined [r8"fT9~] in the third constraint picture. These authors obtained generalized 
paramagnetic susceptibility in the noninteracting regime, and a nonextensive modification 
of the Curie- Weiss law in the context of the mean field model. 

In our current work we make a slight departure from the Refs. [T3p^] . We examine a 
classical arbitrary iV-spin system in a weak background magnetic field without adopting 
the high temperature limit. The thermodynamic quantities in the second and the third 
constraint pictures are evaluated as a perturbative series in the nonextensivity parameter 
(1 — q) by disentangling the g-exponential ( II. 2p . This process of perturbative expansion, 
that may be continued to an arbitrary order, is based on the technique developed in 
[20] and previously used in [2T|22] . In passing from the second constraint picture to 
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that of the third constraint we employ a transformation [22] procedure that allows us to 
express the physical quantities in the later scenario directly in terms of the former. In 
our perturbative expansion for the spins in the background field we retain terms at all 
orders of temperature. Subsequently we study interacting spins in an extended form of 
the mean field model [23] where the field strength coefficient (the proportionality factor) 
is made temperature dependent to accommodate quantum fluctuations among allowed 
configurations [21]. The critical temperature and the Curie- Weiss constant have been 
evaluated in the third constraint framework by retaining terms in the perturbation scheme 
up to the order (1 — q) 2 . In particular, the critical temperature in the nonextensive regime 
increases (decreases) compared to its value given by the Boltzmann-Gibbs statistics for 
the domain q > 1 (q < 1). Our observation qualitatively agrees with the results obtained 
in [T7] where a different definition of temperature for the nonextensive spin system has 
been used. The plan of this article is as follows: Spins in the presence of a weak external 
magnetic field is considered in Sec. [Ill This is followed by the consideration of an extended 
mean field model in Sec. IIII1 Our concluding remarks are given in Sec. IIVI 



II Spins in a weak background field 

The classical Hamiltonian of a system of ./V spins in the presence of a background magnetic 
field H is given by 



N 



E = —fiH cos 9i 



i=l 



(2.1) 

where /i is the magnetic moment of the spins oriented at polar angles (j>i\i — 1, . . . , N) 
with the field. The partition function of the system in the second constraint picture reads 



W) 



n f2n 



,=0 J^ t =0 



N 



1 + (1 - q)/3fiH ^2 cos °i 



i=l 



1=5 N 



J^sin^d^d^. (2.2) 



i=l 



In contrast to the extensive case the available phase space of integration in (12 .2p depends 
on the strength of the magnetic field. In the regime of weak magnetic field 



1 



\l-q\ N 



(2.3) 



the integrand is real and positive in the whole phase space, and the partition function 
may be obtained as follows: 



A' 



Zf {0) = {2ir) N $ 0- N & 1 0), $ = J] [1 + (1 - 



(2.4) 



where the binomial sum &i(x) involving the g-exponentials is given by 



N / \ 

& ^ = J2[J^ n (ex Pg ((iV-2n)x)) A , A = 1 + (1 

n=0 ^ ' 



q)N. 



(2.5) 
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In the extensive q — > 1 limit above sum approaches its well-known classical value &i(x) — > 
2 N sinh^ x. The expression (I2.4p for the partition function exhibits singularities at q — 
1 + - for n — 1, . . . , N. These singularities are similar to the ones observed in [2~Tf2"2] . and 
reflect the fact that the number of degrees of freedom plays a key role in determining the 
allowed range of the nonextensivity parameter q. 

The internal energy in the second constraint picture [3] 

U ?^) = ~Jp (2-6) 



is evaluated by employing the partition function 

uf\(5) =»Hm (6i(/3)r ^ N (iv/r 1 e,0) - e[0)) , (2.7) 

where 

m= (27rY 1 -^ N ^ l - q . (2.8) 

Here and elsewhere the primed functions indicate derivatives with respect to their argu- 
ments. The corresponding magnetization is obtained via the defining relation [13] 



Mf^-^ q Zf\ (2.9) 

On subsequent use of a g-deformed Langevin function L q (x) the magnetization (12. 9p 
assumes the form 

Mf =tiNVlL q 0), (2.10) 

where 

L q {x) = x^ N (©i(x)) 1 ^ (coth g (x; N) - aT 1 ) , coth 9 (x; N) = j^^t - (2.11) 

In the extensive q — > 1 limit, the magnetization (I2.10p reduces to its well-known classical 
value. The magnetization in the second constraint as a perturbative series in (1 —q) reads 

M 9 (2) =/jVI (L0) + (l-q) Mi + (l-q) 2 M 2 + ..X L(x) = coth(x) - x~ l . (2.12) 

The coefficients in the perturbative series up to (1 — q) 2 are given by 

Mi = N 2 L0) In Z0) - N(l - N) /3 2 coth $ cosech 2 /3 + N 2 coth 2 $ 

+N(1 - 2N)/3 coth 2 p-N(l- N)$, 
M 2 = -N 2 (l-2N)/3 + N (I- N)/3 3 -6N (I- N + N 2 )/3 3 coth 2 /3 cosech. 2 /3 

— iV 3 coth p-(N- 5N 2 /3)/3 3 coth 4 p + —(l-N- 12N (3 2 )f3 coth 2 (3 

N 3 - . . . 

+— (5 2 coth 3 (3 + N(3-5N + 2N 2 )(5 A coth 3 (5 cosech 2 /? 

+-0 (( 4 - 13iV + 10Ar2 ) - ( 4 - 6iV + 2iV 2 )/3 2 ) /3 2 coth /3 cosech 2 /3 

7V . TV 3 

+ — (2 - 6iV + 3iV 2 )/3 3 cosech 2 /3 + — (lnZ(/3)) 2 L(/3) 

-iV 3 /3 coth/3 lnZ(/3) L0) + iV 2 (l - JV)/3 cosech 2 /3 lnZ(/3) 
+N 2 /3 2 coth/3 cosech 2 /3 lnZ(/3), (2.13) 
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where Z(x) = 2sinh(x)/x. The perturbative evaluation of the internal energy in the 
second constraint framework is found to follow the standard thermodynamic relation 

U = —M H. (2.14) 

For the sake of brevity we refrain from quoting it explicitly. 

Our task now is to translate the previous results to the third constraint picture. As 
evident from the context of (11.91) the sum of the q- weights plays a seminal role in enacting 
this transformation. The definitions (II. 3p and (11.41) lead to the relation 

c( 2 )(/3) = W) (if))"*. (2.15) 
In the above equation the sum i?(/3) is given by 

f2tf) = [1 " (! " QW^ = (2n) N $ A $-» & 2 0), (2.16) 

i 

where the binomial sum ©2 (x) reads 

6 ^) = E( n )(-l) n (ex P(? ((A-2n)x)) A - (1 - 9) . (2.17) 

n=0 ^ ' 

Employing (12.16P and (12. 4p we may express the sum of g-weights in the second constraint 
picture as 

t {2 \{3) = <R A & 2 0) (6i(/3))-", (2.18) 

and its perturbative evaluation up to terms 0(1 — q) 3 reads: 

c (2) (/5) = k (i + (i - q) Vx + (i - qf r 2 + (i- g) 3 v 3 + ...), (2.19) 

where the perturbative coefficients may be enlisted as 
Vi = N In Zip) — N/3 coth /3 + N, 

P 2 = — /3 2 coth 2 /3 ^(1 — 2N) — 2(1 — N)/3 coth (3 J /3 2 cosech 2 /? 

+ —(^2-2/3 coth/3 + lnZ(/3)J lnZ(/3), 

p 3 = AT 3 /? (1 + 4/3 2 ) coth /3 + A(3 - 5A - 2A 2 ) /3 5 coth 3 /3cosech 2 /3 

7V TV 3 - 

+ — (17 - 33A + 12A 2 )/3 4 cosech 2 /3 - — [3 2 (l + coth 2 (3) 

+ — (17 - 39N + 22A 2 )/3 4 coth 2 (3 cosech 2 /3 

/V 3 - / - ~ „ „ „ \ 

+ — In Z(/3) (3 (3 2 coth 2 (3 + 3(1 - (3 coth /?) In + (In Z(/3)) 2 J 

- N2 ^-~ N \ l - 2 (3 coth (3) (3 2 cosech 2 /3 lnZ(/3). (2.20) 

Following an approach used in [22] we interrelate the dynamical quantities such as the 
internal energy and the magnetization directly from the second to the third constraint 
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picture. The equivalence of the ensemble probabilities f ll.Sp leads to a ready translation 
of the expectation values of an observable O in the two constraint frameworks as 

of (« = Sm- (2 - 21) 

To fruitfully employ this procedure we need to invert the transformation relation (II. 9p 
of the temperature. As a closed form inversion rule is not at hand, we adopt a pertur- 
bative technique [21], and express the auxiliary temperature 0' in terms of the physical 
temperature (3 retaining terms up to second order in (1 — q): 

& = | (1 + (1 - q) g® + (1 - qf h(jj) + ...), (2.22) 

where the rescaled temperature is given by (3 = (3 <! yi~ 1 . The perturbative terms of the 
transformation relation read 

g(f}) = -2N(1 -p coth^) -Nln Z(p), 

— 3iV —i — — s — — — — 

h{fi) = — — (1 - coth 2 (3) - 2N(1 + N)(3 coth (3 cosech 2 /3 - 5N 2 (3 coth (3 

^ N 2 —9 — —9 — —9 N 2 — 

^~ (3 coth 2 (3 + 3iV 2 cosech 2 /? — iV 2 (1 — (3 ) + — (In Z{(3)) 2 
-3N 2 coth^ In Z{fi) + 2iV 2 ^ 2 cosech 2 ^ In Z{fi). (2.23) 

Aided by the inverse transformation series (12.221) we employ (12.211) for obtaining a per- 
turbative evaluation of the magnetization in the third constraint approach to the order 
0((l-q) 2 ): 

Mf ((3) = n (N L(P) + (1 - q) Wit + (1 - qf Wl 2 + • • •) , (2.24) 
where the coefficients OJTi and WI2 read 

j2 i^-q 2 i\ t rS\ t\t2 



SJli = N {2(3 — 1) L(f3) — N (3 In Z(f3) L'(f3) 

N 
~2 



— — N 2 —1. —1 — 1 — — 

Tl 2 = 2N 3 L(f3) + 4iV 3 (3 — (3 — N f3 — (2N + 5N 2 ) (3 coth 2 (3 cosech 2 /3 



+ — (4 + ION - (4 + 9iV 2 ) f3 J (3 coth /3 cosech 2 /3 



+2N(6 + N + N 2 ) /3 4 coth 3 /3 cosech 2 /3 - y (9iV 2 + AN/3 - 8(3 2 )(3 coth 2 /3 
3AT , T ^s AT 2 



(2 + N) /3 3 coth 4 /3 + — (9- 3QN - (4 - QN) f3 2 ^j (3 cosech 2 /? 



2 

+JV/T coth/3 cosech 4 /3. (2.25) 

The internal energy Uq 3 \[3) in the third constraint picture may be read off directly from 
the corresponding magnetization Mq 3 \f3) via the general thermodynamic relation (12.141) . 
We do not quote it explicitly. Magnetic susceptibility is defined as 

xf(P) = (2.26) 
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d 2 = — + — coth (3+ — (N coth 2 (3 - cosech 2 /?) + iV(4 + N + 17N 2 )(3 coth (3 



and we now employ (I2.24p to compute it: 

4 3) (/3) = ((f ) 2 - ^cosech 2 ^ + (1 - q) ^ + (1 - q) 2 § 2 + . . J . (2.27) 

The above perturbative coefficients up to the order 0((1 — q) 2 ) are given below: 

/V 2 — — — 2A^ 2 — 

$i = — 2 -\nZ(f3) + N(l + N + N\nZ((3))cosech 2 (3 + ^ r -2N 2 cosech 2 p 

(3 /3 

— 2 — — 2 — — A^ 2 — 

+N(1 + N)(3 cosech 4 /? + 2iV(l + N)(3 coth 2 /3 cosech 2 /? - coth £ 

—N0 {AN 2 + lnZ(/3)) coth/3cosech 2 /3, 
4iV 3 3N 3 , ^ N 2 

N — — —2 — 

+ y ((6 - + 30N 2 ) cosech 2 /? - (18 + 47N 2 ) coth 2 0) (3 cosech 2 /3 

+4N 2 (N - 1 + (1 + 2A0 coth 2 /3) /f coth/3cosech 2 /3 - N 2 (l - AO/fcosech 2 /? 
+N 3 (7 + 22N) coth /3 cosech 2 /? - N(2 + 3N + 3N 2 ) ~f coth 4 cosech 2 /? 
-iV ((12 + 25N- 19Af 2 )coth/3 + N(3 - 13iV)) /f coth /3 cosech 2 /? 

/V _ _q _ _ 

(40 - 39 - 48N 2 - (32 - 54N + 15N 2 ) coth 2 0) coth (3 cosech 2 ^ 

3 

+iV(2iV(l + N) + (13 - 12N - 7N 2 ) coth 2 /?) cosech 4 /3 

3N 3 - N 3 - - 
+ — 5-ln2:(/3) - ^coth/3 lnZ(0 

-N 2 (l - N) cosech 2 /? In Z(0) - N 2 (10 - JV) /? coth/3 cosech 2 /? lnZ(/3) 
+A^ 2 (5 + 2(1 - 3A0 cosech 2 /? - (19 + ION) coth 2 /3) /f cosech 2 /? \nZ(0) 
+2N 2 ((3 + 7N) coth 2 - 5 AT - 1 + (3 - AOcosech 2 /?) 3 coth cosech 2 /? In 

TV 3 _ 9 /V 3 — — 2 _ _ O 

+ — (In Z{0) - — cosech 2 /? (in + 3A^ 3 (3 coth /3 cosech 2 /? (in Z{(3)) 

2/3 

+ N 3 ( 1 - 3 coth 2 /?) cosech 2 /? (in Z (0) ) 2 . (2.28) 

The specific heat in the third constraint framework is obtained via the corresponding 
internal energy: 

Cf)(/3) = = * ^ M _) - iVcosech 2 /? + (1 - g) ^ + (1 - g) 2 tf 2 + . . . j . 

(2.29) 

The thermodynamic quantities in the weak field limit $ <C 1 follows directly: 



Uf\0 = ~^f^ 2 H 2 , Cj»> = ^ fc/3 2 , (2.30) 



0.98 




1.02 



Figure 1: Dependence of the ratio N eS /N on q for various values of N. 



where the g-dependent effective number of spins N e s reads 



N ( (1 - a)' 

iV eff =^(l-(l-g)(l + iVln2)-L_^L 



(N -2N\n2 + N 2 - (N\n2) 2 ) + . . )j . 



(2.31) 

Substituting the value (12. 8p of the g-dependent scale factor 9T in (I2.3ip . we notice that 
up to the perturbative order (1 — q) 2 it follows N c s > N (N e g < N) for the region q > 1 
(q < 1). This is evident from the Fig. [TJ 

Turning to the magnetization (12.241) in a weak field (3^1 regime, we obtain 



Mf 



N cS fi 2 H 
3kT '■ 



(2.32) 



where the corresponding susceptibility may be viewed as a nonextensive generalization of 
Curie's law: 

(3) _ ^" eff r -Weff^ 

Xq ~ ' 



-cff 
-jT: 



cff 



3k 



(2.33) 



The thermodynamic quantities evaluated above (I2.30H2.33]) show a nonlinear dependence 
on the number of spins N, and the nonextensivity parameter (1 — q). This effect of nonex- 
tensivity manifest in the aforesaid inequality between iV e fj and iV leads to the phenomenon 
of dark magnetism discussed in [T3] . 

The specific heat, the internal energy and the magnetization have also been evaluated 
using the approach discussed in [1], where the sum of the q- weights plays a central role. 
The relation between the ensemble probabilities (ll.Sp enables us to obtain c^ 3 \(3) as a 
perturbative series 



c^(/3) = m (1 + (1 - q) <Px + (1 - q) 2 % 2 + (1 - g) 3 % + •••) 



(2.34) 
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where the coefficients may be listed as 



<Pi = N lnZ(/3) -N (/3 coth/3- 1), 

<£ 2 = ^- (in + AN 2 - 2 /fcosech 2 /? - 2 /3 coth /?) In 

AT 2 / _ _ _9\ _ _ 

+_ (4 _ 4/3 coth /3 - /3 J + N(l + N)f3 coth (3 cosech 2 /? 

AT —9 — 

— (1 + N)p cosech 2 /?, 

/V 3 _ _ _ N 2 —<? —9 _ 

<£ 3 = 5AT 3 -— -/? (/3 + 8coth/3) + — (l-2A^-6A^/3 ) /3 cosech 2 /3 

+^(17 - 12iV + 24A^ 2 )^ 4 cosech 4 ^ - 2N(1 + N + N 2 )ff coth /? cosech 2 ^ 
-JV(3 - iV 2 )/f coth 3 /? cosech 2 /? - iV 2 (l + 3iV)/f coth /3 cosech 4 /? 
-— ~f coth 3 /3 - ^(17 - 23iV - 12N 2 )/? 4 coth 2 p cosech 2 /? + AN 3 In 20) 



6 '4 
— A^ 2 (1 - 2N) /fcosech 2 /? lnZ(/3) - 3iV 3 /? coth/3 lnZ(/3) 

_Q _ _ _ A^ 3 — 9 — — 

-N 2 (3 + 2N)f3 coth/3cosech 2 /3 lnZ(/3) + — /3 coth 2 (3\nZ{(3) 

-N 2 (3 - A0/3 4 coth 2 /3cosech 2 /3 lnZ(/3). (2.35) 

The exponent property (11.71) in conjunction with our perturbative evaluation (I2.34p of the 
sum of the g-weights, now allows us to employ the formulation specified in [4] as a con- 
sistency check on our results. Examining the generalized partition function a differential 
equation involving the internal energy has been established in [I]: 



P dB ~ dB q q dB l-o ' { ] 



where we have used the exponent relation (11. 7p in the last equality. Substituting the 
internal energy that may be readily obtained via the equations (12.141 I2.24p . and the sum 
of the q- weights given in (I2.34p . it may be explicitly verified that the differential equation 
(I2.36P holds order by order in our perturbation theory. This is a nontrivial consistency 
check on our results for physical quantities. It has been noted earlier [22] that, as a 
consequence of the exponent relation (11.71) . evaluations of the thermodynamic quantities 
such as specific heat, say, up to the second order in (1 — q) necessitate computing the sum 
of the q- weights c®(/3) till the third order. That counting has been maintained in (I2.34p . 

Furthermore, the differential equation (I2.36P may be used for a direct extraction of the 
susceptibility in our model. As suggested by (12.241) we postulate the general functional 
dependence of magnetization as 

M® = Hf%. (2.37) 
The relations (I2.14p . (I2.36P and (I2.37P now readily produce the magnetic susceptibility: 
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The derivation following from ( I2.38P precisely agrees with the magnetic susceptibility 
obtained earlier in (I2.27p . This confirms the validity of our perturbative procedure. 



Ill Mean field model: temperature dependent 
effective field coefficient 

Physical reasoning tells us that effects of nonextensivity is likely to be pronounced for 
systems embodying long-range interactions between the constituents. For an interacting 
spin system a good first approximation is provided by the mean field model where the 
long-range component of the interaction between the spins is taken into account via a 
homogeneous magnetic field (H m ) that is assumed to be directly proportional to the 
magnetization per spin. The Hamiltonian of the system reads 



N 



E = -fi(H + H m ) ^ cos 9i 



(3.1) 



i=i 



where H m is the dynamical field resulting from the long range interactions between the 
spins as envisaged in the mean field model. The generalized partition function in the 
third constraint reads 



Zf 03) 



exp„ [e) 



2tt 



J<f>i=0 



N 



1 + (1 



q)P ^2 cos 0i 

i=i 



— N 

Y\ sin 6i d9i dfa 

i=l 



(3.2) 



where e 



(3)-, and the scaled dimensionless variable (3 is given by 







c V)(p) + (l-q)pU, 



(3)' 



(3-3) 



In the regime |1 — q\ < 1 we assume that the integrand is real and positive in the entire 
phase space, and obtain the generalized partition function as follows: 



Zf (ft) = (270* $ exp q (e)p- N S^/3). 



(3.4) 



The primary definition of a thermodynamic observable fll.5p leads to the following 
integral form of the internal energy in the third constraint picture: 



uf\p) 



f,(H + H m ) (Z^ifS))- 1 (exp 9 (e)) 



/7T p27V N 

_ i= o J4, i= o i=l 



N 



l + (l-q)pY, cos0 i 



i=l 



JJsinMW;. (3.5) 



i=l 



Implementing the above phase space integration we obtain the internal energy as 



uf\p) 



2tt 
J 



N 



<S>fi(H + H m ) (n/3- 1 ©!(/?) -e[(p) 



(exp g (e)) 1 



K3) 



(3.6) 
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On substituting the generalized partition function ( 13. 4 p the implicit equation (13. 6 p may 
be recast as follows: 



fx(H + H m ) IN &M 
1 + (l-q)e\p <S X 0) 



(3.7) 



We expand the term in the parenthesis up to second order in the nonextensivity parameter 
(1 — q) and fourth order in the dynamical variable Pfi(H+H m ). The nonzero contributions 
in the resulting expansion read 



N(3fi 2 (H + H n 



3 c( 3 )(/3) 

where the perturbative coefficients are given by 



_lk f(3fx(H + H v 
1 15 



c( 3 )(/3) 



(3.f 



Ui = 1 - (1 -q) (1 + 2e) + (1 -q) 2 e(2 + 3e) + 
U 2 = l-(l-g) {6-10N + ie) 

+ (1 - qf (11 - 2hN + 2Ae - 40e + lONe 2 ) + .... 



(3.9) 



The phase space integral corresponding to the sum of the g-weights may be read off via 
(USD and flUD: 



c (3) (/3) 



2tt 



i= J<t>i=0 



N 



1 + (1 -q)P^2cos0i 



i=l 



i 

1-s N 



J^sin^d^d^. (3.10) 



i=i 



Performing the above integrations and subsequently substituting the generalized partition 
function (I3.4p above sum of the g-weights assumes the form 



-(3) 



{P) = m a (5 {q ~ 1)N e 2 (f3) (e^))-". 



(3.11) 



As done before in the instance of the internal energy in (13. 8p the rhs of the expression 
(13. lip is expanded perturbatively up to second order in the nonextensivity parameter 
(1 — q) and fourth order in the variable f3fi(H + H m ): 



c (3) (/3) 



m(l + (1 - q)N\n2 + (1 - q) 2 y (l + N + N(\n2) 2 ) 



+Pi(e, 



V c(3)(/3) 



+ P2(e) 



Pfj,(H + H n 
c( 3 )(/3) 



+ 



(3.12) 



where the perturbative expansion of the coefficients of monomials of the field variable 
read 



N 



6 



((1 - q) + (1 - q) 2 (1 - Win 2 + 2 e) + . . .) 



N 



P2(e) = ^[(l-q)-(l-q) 2 [Q 



15N 



AHn2 + 4£ 



(3.13) 
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Towards obtaining the magnetization we first explicitly obtain the quantities Uq (/?) 
and c^ 3 \(3) by solving the pair of simultaneous implicit equations f 13 . 8 j) and f)3.12p . We 
can systematically obtain their solutions in an order by order perturbation theory where 
we retain terms up to (1 — q) 2 in the nonextensivity parameter and ((3 fi(H + H m )) A in the 
field variable. The relevant expression for the internal energy reads 

rr(3)fm _ N ^\H + H m f N gV%g + E m f u 
q {P) 3 m + 45 W - + K^V 

where the g-dependent numerical factors T and H are given by 

T = 1- (l-g)(l + JVln2) 

-(1 -g) 2 y (l + JV-2 ln2-iV(ln2) 2 ) + ..., (3.15) 

S = l-(l-g) (6+™+3N\n2 

+ (1 - qf (\l - ^(A^ + 1) + (36 + + 9iVln2) )+.... (3.16) 

The perturbative expansion also yields the sum of g-weights as 

,(3) 



m(^l + (1 - g)AHn2 + (1 - qf y (l + N + iV(ln2) 2 ) 

+n ^m^\ n2 (m±Mj + ..^ ( 3, 7 ) 



where the coefficients may be listed as follows: 

n, = _^((l- g ) + (l- g ) 2 (l-iVln2) + ...), 

n 2 = ^^l-q)-(l-qf(^ + ^ + 3N\n2 ) )+..)). (3.18) 

As the internal energy and the magnetization are related via the standard thermodynamic 
expression (12.141) we may now readily obtain a perturbative expansion for the magnetiza- 
tion by employing the corresponding series f !3.14p for the internal energy: 

=t*(J^ + H m ) r — {^{H + H m )f e) . (3.19) 

The homogeneous magnetic field if m is taken to be proportional to the magnetic moment 
per spin. To account for the quantum effect of spin fluctuations the proportionality factor, 
known as the effective field coefficient (A(T)), is considered to be temperature dependent 
[2~4"f2"5] . Following [23] we here study a mean field spin model with the field coefficient 
A(T) being linearly dependent on temperature. The above discussion leads to the field 
variable 

M (3) 

H m = \(T)mf ) rn^ = ^, \{T)=i + CT ) (3.20) 

12 



where the coefficients £ and ( characterize the long range spin interaction. Substituting 
(I3.20p in (I3.19P and considering the vanishing limit of the the external field H = 0, the 
magnetization reads 

-p 9 (w + f )' 92 (? ) s «') 2 + 1 - ¥ i = °< (3 - 2i > 

where the critical temperature is 

To visualize the variation of the critical temperature with respect to the nonextensivity 
parameter (1 — q), we, in FigF2J plot the ratio k defined as 

From Figf2] it may be inferred that the critical temperature increases (decreases) com- 
pared to its standard Boltzmann-Gibbs value in the regime q > 1 (q < 1). This result 
is qualitatively similar to the observation in [17] where these authors have adopted an 
alternate definition of the temperature for the nonextensive spin system. The mean field 
model studied by these authors is slightly different from the one considered here in that we 
assume the dynamical field strengh to be temperature dependent in order to accommodate 
the quantum fluctuations in spin configurations. 

Equation (13.211) suggests that above the critical temperature T > T c ■ , the only real 
solution for the magnetization in null external field condition is given by 



in 



fk=o = (3.24) 



that corresponds to the paramagnetic phase. Using the standard definition (I2.26P in 
conjunction with the relations (I3.19P and (I3.24p we now obtain the magnetic susceptibility 
in the paramagnetic regime: 

(3) _ C r - N ^ T 

x q - r _ r ( 3 )> c -^kme- (3 - 25) 

It is evident from (I3.2ip that below the critical temperature T < the magnetization 
in the null external field limit is given by two stable real values given by 



Ik p3 / /■ \ ~3/2 j, / 

<VM» = ±A*V3g e-^l + |Tj (3.26) 



This corresponds to the ferromagnetic transition as observed in the generalized mean field 
approach in the nonextensive framework. In the regime T — > T^f 1 ~ the susceptibility in 
the null external field limit may be obtained via (I3.19P and (I3.26p : 



4" = (3-27) 



13 



0.98 




1.02 



Figure 2: Variation of the quantity k with q for various values of N. 



As characteristic of the mean field approach, the divergence of the susceptibility follows 
critical exponent law 



4 3) ~ 7^7^)? ( 3 - 28 ) 



1 

|T-T C V 

where the critical temperature depends on the nonextensivity parameter (1 — q). Of 
course in the extensive q — > 1 limit, the standard physical quantities corresponding to the 
Boltzmann-Gibbs statistics are recovered. We also quote the results obtained in the usual 
mean field model, where the dynamical proportionality factor A is regarded as independent 
of temperature i.e. the linear coefficient in f)3.20p assumes the value ( = 0. In this limit 
the critical temperature and the Curie- Weiss constant in the nonextensive scenario may 
be arrived at via ( 13 . 22 j) and (13. 25ft : 

^ £At 2 r Nu 2 T , 

T ^-jk- £ =W' < 3 29) 



IV Remarks 

Our main focus in the present work has been to study a system of spins with long range 
interactions approximated by a mean field model governed by the nonextensive Tsallis 
statistics. To incorporate the quantum spin fluctuations the mean field model investi- 
gated here is assumed to have the dynamical field strength coefficient depending linearly 
on temperature. The nonextensivity is implemented by using a perturbative technique, 
where the implicit simultaneous equations involving the internal energy and the sum of 
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g- weights were solved explicitly as series expansions up to the order (1 — q) 2 in the nonex- 
tensivity parameter. The perturbation method developed here may be continued to an 
arbitrary order in the parameter (1 — q). The signature of the nonextensivity is evident as 
the critical temperature is found to depend on the number of spins N and the deformation 
variable (1 — q). Compared to its standard Boltzmann-Gibbs value the critical tempera- 
ture increases (decreases) for the domain q > 1 (q < 1). The extended Curie- Weiss law 
characterizing the susceptibility in the regions above and below the critical temperature 
has been generalized to the nonextensive case. Analogous to the critical temperature 
the Curie- Weiss constant also embodies the effects of nonextensivity. Another interest- 
ing feature reflecting nonextensivity emerges for noninteracting spins in the presence of a 
background field, where we observe the presence of dark magnetism. This supports the 
results [T37T4"] obtained earlier. Parallel to the observation in [T4] our analysis indicates 
that the effective number of spins N e ff > N (N c g < N) for the regime q > 1 (q < 1). 

The mean field technique used in the current work may be fruitfully employed in 
calculating the magnetic properties of systems which form clusters [25]. A cluster is 
typically a macroscopic region consisting of a large number of spins whose interactions are 
described via locally homogeneous mean fields in the domain of each cluster. Then each 
cluster is approximated by a single effective spin forming an ensemble of varying spins 
whose magnetic properties may be studied by adopting a suitable model. The results 
obtained here may facilitate studying such models in the context of Tsallis statistics. 
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